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ABSTRACT 

We present a new method aiming to simplify the cosmological analysis of X-ray clus- 
ter surveys. It is based on purely instrumental observable quantities, considered in a 
two-dimensional X-ray colour-magnitude diagram (hardness ratio versus count-rate). 
The basic principle is that, even in rather shallow surveys, substantial information on 
cluster redshift and temperature is present in the raw X-ray data and can be statisti- 
cally extracted; in parallel, such diagrams can be readily predicted from an ab initio 
cosmological modeling. We illustrate the methodology for the case of a 100 deg2 XMM 
survey having a sensitivity of ~ 10 ergs/s/cm 2 and fit at the same time, the survey 
selection function, the cluster evolutionary scaling-relations and the cosmology; our 
sole assumption - driven by the limited size of the sample considered in the case-study 
- is that the local cluster scaling relations are known. We devote special care to the 
realistic modeling of the count-rate measurement uncertainties and evaluate the poten- 
tial of the method via a Fisher analysis. In the absence of individual cluster redshifts, 
the CR-HR method appears to be much more efficient than the traditional approach 
based on cluster counts (i.e. dn/dz, requiring redshifts). In the case where redshifts are 
available, our method performs similarly as the traditional mass function (dn/dM /dz) 
for the purely cosmological parameters, but better constrains parameters defining the 
cluster scaling relations and their evolution. A further practical advantage of the CR- 
HR method is its simplicity : this fully top-down approach totally bypasses the tedious 
steps consisting in deriving cluster masses from X-ray temperature measurements. 

Key words: Galaxies: clusters: general - Cosmology: observations - X-rays: galaxies: 
clusters - Methods: observational 



1 INTRODUCTION 

Discriminating between different cosmic scenarios requires 
precision cosmological studies relying on well-controlled 
observables. In parallel to the Cosmological Microwave 
Background (CMB), Baryonic Acoustic Oscillations (BAO), 
Type la Supernovae and Weak-Lensing analyses, galaxy 
clusters, as the most massive bound entities in the uni- 
verse, are expect ed to provide independen t comp l emen- 
tary constraints 
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particular, they appear to be quite sensitive to the prop- 
erties of the dark ene rgy JHaiman Mohr, fc Holdeil l200ll ; 
iBattve fc Welled 120031 ; iPierre et alj l201ll ) Cluster cosmo- 
logical studies are usually based on the cluster number 
counts as a function of redshift and mass. This quan- 
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tity can be easily inferred from the halo model for- 
malism and is confirmed by the most recent N-body 
simulations (IPress fc Schechterl 1 1974 ILacev fc Cold Il993; 
Sheth fc TormenI Il999l; Ijenkins et al.1 l200ll ; ISpringel et al.1 
20051 ) ; see lCoorav fc Shethll2002l foTareview. Theoretically, 
its high sensitivity to the initial density fluctuation power- 
spectrum as well as its time-evolution make it a powerful 
probe of structure formation. Its dependence on geometrical 
effects (surveyed volume) further strengthens the constraints 
on key cosmological quantities such as the matter content 
in the Universe (f2 m ) and the dark energy equation of state. 
From the observer's point of view, however, cluster masses 
are not quantities easily measurable (contrary to cluster red- 
shifts), a fact that often leads to question the actual use of 
clusters as cosmological probes. 

Galaxy clusters can be studied in a variety of ways, in 
particular through their X-ray emission. The gas trapped 
in the deep cluster potential is heated up to X-ray emit- 
ting temperatures. Free-free emission is the dominant mech- 
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anism from the hot plasma having a heavy element abun- 
dance of ~ 0.3 Zq; at low temperatures (<2 keV) a sig- 
nifica nt fraction of the energy is emitted via recombination 
lines l|Sarazinll 19881 ). Because extended X-ray sources at high 
galactic latitude almost unambiguously point toward clus- 
ter potential wells, hence minimising projection effects, X- 
ray surveys have long been considered as the ideal way of 
constructing cosmological cluster samples. 

The Einst ein Observatory Extended Medium Sensi- 
tivity Survey (iGioia et al.l [1990) provided the first flux- 
limited, X-ray selected sample of galaxy clusters, allow- 
ing pioneering clust er count analyses. It was then fol- 
lowed by REFLEX (TBohringer et al.l l200ll ) and NORAS 
jBohringer et all 12000! ) based on the ROSAT all-sky sur- 
vev (lTruemperlll993h as well as a by a numb er of cluster 
searches in the deep ROSAT archival po i nting s (Schar fet all 
1 19971 ; iRosatill 19981 ; Ivikhlinin et~allll998l ; iRomer et al.ll2000l ). 
Preliminary cosmological constraints resulted from these 
studies involving not o nly cluster counts, bu t also their 3- 
D spatial distribution (jSchuecker et al.ll2003h : at this stage 
however, the main observable quantity that was dealt with 
was the c luster luminosity function, rather than t he mass 
function l|Borgani et all 1200 ll ; iMantz et ail I2008T ). Paral- 
lel analyses invoking the distribution of temperatures in 
R OS A T and ASCA clusters also pro vided cosmological con- 
straints, however somewhat debated dHenrvll 19971; lEke et al.l 
199Sl;IViana fc Liddle|[i999l; iPierpaoli. Scott, fc Whitell2o5lT 
Pierpaoli et~ail 120031 ; iHenrvl |2004 iHenrv et aill2009h . but 



consistent with findings from the X-ray luminosity func- 
tion. One of the major shortcomings of these studies rapidly 
turned out to be lack of reliable mass- observable relations, 
ideally in the form of scaling relations, and how these would 
evolve as a function cosmic time. 

With the advent of XMM-Newton and Chandra, previ- 
ous samples underwent deep observations and, in parallel, 
the interest in X-ray surveys for cosmological analyses in- 
creased. A very significant amount of observing time was 
devoted to the determination of the cluster scaling rela- 
tions, for samples a priori th ought to be repre sentative of 
some cluster population (see iPratt et alj (I2009T) and refer- 
ences therein). In particular, the 400d survey ( Burenin et al.l 
120071 ) provided a sample of some 90 clusters selected in the 
RASS and followed-up by deep Chandra observations, allow- 
ing precise mass measurements based on hi gh-quality X-ray 
data. T he cosmological analysis presented in lVikhlinin et all 
(2009b) relies on the mass and redshift distribution of clus- 
ters. lMantz et all (|2010al ) used more than 200 X-ray selected 
clusters, some of them having a deep follow-up, in an anal- 
ysis that combines the gas mass fraction in clusters with 
their abundance per mass and redshift bin. Simultaneously, 
new cluster samples have been assemble d either from dedi- 
cated XMM surveys (|Pierre et al.ll2004) o r from the XMM 
and Chandra ar chival data: e . g. ChaMP dBarkhouse et alj 
I2006D and XCS (|Romer et all 1200 it iMehrtens et alll201ll ). 
Corresponding cosm ological analyses are still in progress 
l|Sahlen et all 12009). but one of the major outcome was to 
realise that selection effects can be as critical as the proper 
knowl edge of the cluster scaling relations. iPacaud et all 
|2007l) have shown that, unless a high flux limit is as- 
sumed, X-ray cluster samples are best characterised by a 
two-dimensional selection function (analogous to a surface 
brightness limit). Further, because of the steepness of the 



cluster mass function and of the (currently poorly deter- 
mined) dispersion in the scaling laws, these relations ap- 
pear to be always biased toward the most luminous objects 
with respect to the mean, unless a thorough treatment of 
the selection function is introduced; this becomes especially 
challenging as redshift increases. Ideally, for a given clus- 
ter sample, one would need to simultaneously model (i) the 
selection effects (ii) the scaling relations and (iii) the cos- 
mology. 

From the observer's point of view, the bottle neck in 
building large cosmological samples is the time-consuming 
optical follow-up to obtain spectroscopic redshifts for each 
cluster. However, redshift (and mass) information is already 
encoded in the X-ray spectra of the clusters. In principle, 
it is possible to make use of this information in a statis- 
tical way, even i n the low-count regime, thanks to a ded i- 
cated formalism l|Llovd-Davies et alllioill ; lYu et aill201ll ). 
The goal of the present article is to investigate such a new 
approach to the cosmological analysis of large samples of 
X-ray clusters: the CR-HR method. Conversely to methods 
requiring redshift information for each cluster and inferring 
the cluster mass distribution through various X-ray proxies, 
we handle X-ray instrumental observables only namely, the 
count-rates in several energy bands. In the scientific analy- 
sis, we self-consistently model the cosmology, scaling laws, 
selection effects and the instrumental responses to predict 
count-rate distributions that can be directly compared to 
the purely observational data. 

The structure of this paper is as follows. We begin by 
presenting the motivations of the CR-HR method and de- 
scribe its principle. We then give the key ingredients involved 
in the construction of the CR-HR diagram that we illustrate 
for a shallow XMM survey (Sect. [3}- Then, we explain the 
modeling of measurement errors and their inclusion in the 
analysis (Sect. [4]). Next, we describe the formalism adopted 
in our Fisher analysis used to evaluate the CR-HR method; 
we present the expected constraints for a set of selected pa- 
rameters (Sect[5]). We discuss and summarize our results in 
Sect. H 

Throughout the paper, we assume a flat (Qk = 0) 
ACDM co smology with param eters given by WMAP-5 best 
fit values (iDunklev et al.l l2009). 



2 THE CR-HR METHOD 

In this work, we consider a shallow X-ray survey and as- 
sume that a robust procedure allows the construction of 
well-defined samples of clusters of galaxies. By definition, the 
survey selection function is based on X-ray observable crite- 
ria only. The survey is supposed to be shallow in the sense 
that a few hundred photons, at most, are collected for each 
cluster and, thus, may enable a mean temperature estimate 
of the intra-cluster medium but no radial temperature pro- 
files. This is the case f or most of the current analyses to date 
llPacaud et ail (l2006l). ISahlen et all (|2009l ). lBarkhouse et all 
(2006). iBurenin et all 720071 )). Given a cluster sample, we 
populate a 2D observable parameter space defined by the 
measured X-ray count-rate (CR) and hardness ratio (HR). 
CR and HR, which contain (partially degenerate) informa- 
tion on the temperature and the redshift of the clusters, are 
defined for adequately chosen X-ray bands. Specifically, we 
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construct a CR-HR diagram from the selected cluster sam- 
ple, building a two-dimensional density, which behaves like 
an X-ray color-magnitude diagram. In the case where op- 
tical (photometric or spectroscopic) redshifts are available 
for each cluster, we divide the sample in redshift bins and 
associate a CR-HR diagram to each of these slices leading 
to a three-dimensional z-CR-HR diagram. 

Such a diagram can, in turn, be obtained using an ab 
initio formalism: (i) setting a cosmological model, we com- 
pute the number of clusters as a function of mass and red- 
shift; (ii) each cluster is ascribed an X-ray temperature and 
luminosity as well as a physical characteristic size according 
to empirical scaling laws; (iii) these quantities are subse- 
quently converted into CR, HR for a given X-ray survey 
and into an apparent size; (iv) finally, only clusters passing 
the X-ray selection function are retained, enabling the con- 
struction of theoretical CR-HR diagrams. We compute such 
(z)-CR-HR diagrams for a wide range of cosmological mod- 
els and possible cluster evolutionary scenarios in order to 
determine which one is the most likely, by comparing with 
the observed diagram. 

Count-rates are purely instrumental quantities, thus in 
order to test the ability of the CR-HR method to constrain 
both the cosmology and cluster evolutionary physics, we 
need to explicit the calculations for a given X-ray survey 
instr ument such as XMM-Newton or eRosita (Prede hl et al.l 
2010) for instance. In the present article, we assume a 
100 sq. deg. XMM survey performed with the EPIC instru- 
ments with a sensitivity of ~ 10~ 14 ergs/s/cm 2 in the [0.5-2] 
keV band for cluster-type sources, i.e. consisting of 10 ks ex- 
posures. The CR-HR method is evaluated by a Fisher anal- 
ysis and its efficiency compared to the traditional method 
relying on the redshift-mass distribution of clusters. Special 
care is given to modeling of measurement errors in the Fisher 
analysis. 



3 INGREDIENTS ENTERING THE CR-HR 
METHOD 

3.1 Modeling the CR-HR distribution of clusters 

This section describes the steps entering the computation of 
CR-HR diagrams. We illustrate our method under realistic 
conditions and qualitatively show its sensitivity to relevant 
parameters entering the model. Measurement errors are con- 
sidered in Sect. [J A schematic view of the method is given 
on Fig. [T] 



3.1.1 Cosmological mass function 

We start from a scale-inv ariant primordial spectr um with 
slope n 3 — 0.961 and the lEisenstein fc HvJ (1 19981 ) transfer 
function to obtain the z = power-spectrum P(k,z = 0), 
which is subsequently normalized by erg. The linear power 
spectrum P(k, z) is evaluated using the redshift-dependent 
growth factor computed by numerical integration of the par- 
tial differential equation. We then compute the rms vari- 
ance <j(M, z) of the field smoothed at a comoving scale 
R = (3Af/47rp m ) 1,/3 and inject it into the following func- 
tional form describing the differential comoving density of 



haloes per mass interval dM about M at redshift z: 

dM ~ ;{a> M dM U 

where p m is the mean matter density at redshift z. We cal- 
culate the mass function in terms of A'^oob, the mass within 
a radius i?200b, inside which the mean mass density is 200 
times the ma tter density in the U niverse 

We use iTinker et all (|2008h fit for to obtain f(a) 
for the corresponding mass definition and then com- 
pute the sky-projected, redshift-dependent mass function 
dn/dfl/dMwob/dz. The equation of state of dark energy is 
parametrized through a single parameter wq = P/p, whose 
value in the case of a cosmological constant is —1. 

We further transform the mass function in terms of 
Afeooc, defined relatively to the critical density of the Uni- 
verse. This conversion (M2oob to M2ooc) is motivated by the 
fact that M200C is the mass definition entering our scaling- 
law formulae (see Section [3l) and is perfo rmed using the fit- 
ting formula from lHu fc Kravtsovl 1I2003I'); for this purpose, 
we as sumed a NFW mass profile dNavarro, Frenk. fc W hite 
Il997h and a concentration model from Bullock et al. I <|200ll) . 



3.1.2 Cluster X-ray emission: scaling laws and brightness 
profiles 

The X-ray emissivity of clusters basically depends on 
three quantities: the redshift z, the cluster X-ray tem- 
perature T integrated over the whole cluster extent, 
and its total bolometric luminosity Lx (along with 
some dependence on the metallicity of the ICM). Scal- 
ing relations between cluster masses and these quanti- 
ties have been extensive l y studied in the local Universe 
jArnaud fc Evrard 19991: Arnaud. Pointecouteau. fc PrattJ 
120051 : IVikhlinin et all 120061: iPratt et alj|2009| 1 down to the 
low-mass end (|Sun et al 1 l2009l V As usual, we model the clus- 
ter scaling relations by power-laws. Given that the physical 
processes determining the evolution of these relations are 
still a matter of debate, we param etrize the evolution by 
the factor (1 + z) 1 (e.g. | Voitj [20051) . Our mass-observable 
relations read: 



M200C 



lO^/i-iJVlo 



n f T \ «MT 
_ "J^Q^MT / 1 



UkevJ 



£(z)- 1 (l + *) 7 -- MT (2) 



10 44 ergs/s 



= 10 



,4keV/ 



e(z)(i + z y 



(3) 



Further, the intrinsic scatter in those relations is an impor- 
tant ingredient for the modeling the c luster population (e.g. 
IStanek etT al. 2006 P acaud et al]|2007h . We introduce a inT \ M 
and cr ini | T , the scatter in T at fixed M200C (respectively in 
Lx at fixed T) and assume they are independent of redshift, 
mass and temperature. 

Finally, we a ssume a surface brightness profile given 
by a /?-model l|Cavaliere fc Fusco-Femianol [l976) with 
j3 = 2/3 and a varying core radius r c . The scaling 
of r c with other cluster quantities is complex and de- 
pends o n the details of the intra-cluster medium physics 
(see e.g. Sanderson fc Ponmanl [20031 . lota fc Mitsudall2004 
lAlshino et al.ll2010h but it can reasonably assumed that r c 
scales with the size of the dark matter halo. We thus take 
a x c fi = r c /i?5ooc parameter, constant at all redshift and 
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Figure 1. Schematic view of the CR-HR method (right part of the block diagram), illustrating the top-down approach used to link X-ray 
observabies to a cosmological model. Left part of the block diagram shows the more traditional method based on individual cluster mass 
measurements using, e.g. a temperature proxy (Tx) and a mass-temperature scaling relation (Tx — M2ooc)- The C symbol indicates the 
comparison between model and data (based on, e.g., a minimization of the likelihood function): the CR-HR method compares directly 
X-ray observabies. 



masses (with iisooc being denned as the radius enclosing a 
mean density of 500 times the critical density of the Uni- 
verse). 



3.1.3 An instrumental model for XMM observations 

Most of the cluster detection algorithms in the X-ray wave- 
band are based on a two-step procedure: source detection is 
run on a filtered image, followed by fitting a cluster emis- 
sion model on the raw photon ima ge, accounting for the 
Poissonian nature of the signal (e.g. iBohringer et alj 200 ll: 
Burenin et al.ll2007l ; iPacaud et al.1 120071 : iLlovd-Davies et al.l 
201ll ). The efficiency of such an algorithm, in terms of com- 



pleteness and purity, is evaluated by extensive image simula- 
tions. This finally enables the determination of cluster selec- 
tion functions based exclusively on X-ray criteria, which are, 
in general, more com plex than a simple flux limit. Following 
IPacaud et all (|2006t ). we use a two-dimensional parametrisa- 
tion involving the count- rate in the [0.5-2] keV band and the 
apparent core radius. Figur e [2] shows our adopt ed selection 
function (the CI selection. IPacaud et al.l [20061 ) which cor- 
responds to an uncontaminated cluster sample, for XMM 
exposure times of the order 10 ks. 

Each cluster is characterized by a redshift z, a tempera- 



ture T and a bolometric luminosity Lx ■ Count-rates are de- 
rived from physical fluxes for a given spectral emission model 
and using the proper instrumental responses. I n this work, 
we a ssume a thermal plasma model (APEC , ISmith et al.1 
l200ll ) having a metal abundance of 0.3 Zq along with a 
galactic absorption corre sponding to Nh = 3 x 10 20 cm -2 
IIGrevesse fc Sauvallll998h . Response matrices for the three 
EPIC detector s onboard XMM (M OS1, MOS2. lTurner et all 
120011 : and PN. IStriider et al.ll200ll ) and THIN filter are used 
to produce an observed spectrum (number of counts col- 
lected by second in each energy channel) which is in turn 
integrated over specific energy bands to yield the desired 
countrates. 

For the purpose of our analysis, we define three work- 
ing energy bands: [0.5-2] keV (band "tot", which is also the 
detection band), [1-2] keV (band "1") and [0.5-1] keV (band 
"2" ) . The choice of these bands is documented in Appendix 
B. In all what follows, we assume that we measure total 
countrates, i.e. over the full cluster extent (see discussion in 
App.P 

We define the cluster hardness ratio by HR = 
CR1/CR2 where CRi and CR2 are count-rates measured in 
[1-2] and [0.5-1] keV respectively. We neglect possible spa- 
tial variations of the hardness ratio across the cluster X- 
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Figure 2. Selection function used throughout this analysis, ob- 
tained from realisti c simulations of XMM cluster observations 
llPacaud et al,l [20061) . The detection probability is expressed as 
a function of two observable quantities: the total count-rate col- 
lected by the three detectors, and the core radius of the input 
/3-model (/3 = 2/3). 
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Figure 3. Lines of iso-hardness ratio CR1/CR2 in the plasma 
temperature-redshift plane for the XMM EPIC instrument. CRi 
is the count-rate in [1-2] keV and CR 2 in [0.5-1] keV. An APEC 
plasma model with abundances 0.3 Zq and a galactic absorption 
Njj = 3.10 20 cm" 2 is used. Changes in abundance value from 0.1 
to 0.6 Zq are indicated by the dotted coloured lines. 



ray extent. In practice, for the type of surveys considered 
here, both the faintness and the small extent of the objects 
(compared to the instrumental PSF) prevent from resolving 
such detailed structure. We thus treat the cluster emission 
as equivalent to that of a single-temperature plasma and 
this is consistent with the fact that reference scaling rela- 
tions have been computed by fitting a single plasma model 
to various cluster spectra. For a given spectral model de- 
pending on parameters (z, T, Lx), the hardness ratio does 
not depend on luminosity. Figure [3] shows the redshift and 
temperature dependence of the hardness ratio values : a pure 
bremsstrahlung spectrum would exhibit a degeneracy of the 
T/(l + z) type; however the presence of metallic lines en- 
tering the energy bands at different redshifts induces more 
subtle effects, especially at low temperatures where they are 
prominent. 



3.1.4 Summary of the model parametrization 

Table Q] summarizes the main parameters used in our anal- 
ysis. Beside the WMAP5 cosmological model, parameters 
governing the cluster M-T and L-T scaling relations are de- 
fined by_ejp 1 sJ5]_amiJ31_r^ T rela tion was taken 
from lArnaud. Pointecouteau. fc PrattJ ([2005) usi ng their re- 
lation for hot clusters and <5 = 200. Following I Alshino et al.l 
we set x c ,o =0.1, the ratio between the core radius 
of the X-ray / ?-model and Rsooc- The local Lx — T relation 
is taken from iPratt et al.l {2009) using their L\ — T\ rela- 
tion for "Non Cool Core clusters". We justify this choice 
by finding that our fiducial model along with the selection 
function of Fig. [5] yields 5.7 clusters per sq. deg., consistent 
with the ob served density of c lusters in the 10 ks XMM- 
LSS survey jPacaud et al.l l2007). Choosing their relation for 
"All" clusters would lead to a higher density of 10.6 clus- 
ters per sq. deg, indicating an incompatibility between this 
relation and the XMM-LSS selection function, possibly orig- 
inating from the different fractions of cool core clusters in 
the samples under study. A typical cluster at z = 0.4 with 
M 20 



1O 14 /i _1 M has a temperature T = 2.2 keV, a to- 



tal bolometric luminosity Lx = 0.55 x 10 44 ergs/s, and radii 
i?500c = 0.6 Mpc and r c = 11 arcsec; with a 10 ks XMM 
exposure, we collect ~ 500 photons for this object. 



3.2 Illustrative examples 

We show on Fig. [4] the CR-HR distribution computed for 
our fiducial set of parameters and illustrate the effect of a 
parameter change on this diagram. The most obvious ef- 
fect of modifying one of the parameters is a variation in 
the total number of observed clusters. This is particularly 
striking for Sl m and as which strongly impact the ampli- 
tude of the distribution and thus are relatively well con- 
st rained by the total number count alone (as pointe d out 
in lHaiman. Mohr. fc Holded l200ll . ISahlen et all 120091 '). Pa- 
rameters governing scaling laws (7^ mt, 7z,lt) enter at the 
cluster selection stage: increasing e.g. 7 z ,lt increases the 
luminosity and thus the detectability of clusters at higher 
redshifts. Beyond this first-order overall change in ampli- 
tude, the shape of the distribution is affected in various 
ways when the model parameters are varied. If one is able to 
detect these changes within the measurement uncertainties 
and systematic errors, degeneracies between parameters can 
be broken. For instance, a change of 20% in the value of Q m 
uniformly changes the amplitude of the CR-HR distribution 
while a +1 modification in 7z,lt also shifts the center of the 
distribution towards lower HR and lower CR. 



3.3 Generalizing: adding redshift information 

We also consider the case where the survey benefits from 
an optical spectroscopic follow-up, providing cluster red- 
shifts. To model this case, we define thin redshift slices 
and repeat the procedure described above in each of these 
slices to derive the corresponding tridimensional quantity 
dn/dz/dCR/dHR. The resulting distributions are illustrated 
on Fig. [5] for given redshift ranges. Such diagrams almost 
fully characterize the whole cluster population using purely 
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Parameter 


Fiducial value 


Description 


Rcf. 


Prior 




0.249 




1 


No 




1 - o m 


(Flat Universe) 


- 


- 


n b 


0.043 




1 


0.003 


(Jg 


0.787 




1 


No 


IV () 


-1 






No 




0.961 




1 


0.014 


ft 


0.72 




1 


0.026 


CI \] •p 


1.49 


M — T power-law index 


2 


0.17 




0.46 


A^f — T logarithmic normalization 


2 


0.023 


T^MT 





Af — T evolution index 




No 


a lnT\AI 


0.1 


M — T constant logarithmic dispersion 


2 


0.064 


qlt 


2.89 


L — T power-law index 


3 


0.21 


Clt 


0.40 


L — T logarithmic normalization 


3 


0.026 


7z,LT 





L — T evolution index 




No 


°"inL|T 


0.267 


L — T constant logarithmic dispersion 


3 


0.058 


^c,0 


0.1 


/3-model core radius scaling wrt. i?500c 


4 


No 



Table 1. Fiducial para meters used in this study. Last column shows the sta ndar d priors used in the Fisher analysis (see sect. (5J. 
References: (D: lDunklev et al.ll2009l , (21: lArnaud. Pointecouteau, fc Pratj|2005l . (3): |Pratt et al.ll2009l , (4): lAlshino et al.|[20icl 



observable (instrumental) quantities: the redshift distribu- 
tion, the evolution of the CR-HR distribution as a function 
of redshift and how these quantities are related in the sam- 
ple. In fact, the dn/dz/dCH/dRK distribution is analogous 
to the dn/dz/dLx /dT distribution but can be readily ob- 
tained from the available data without any assumption on 
scaling laws and on cosmological parameters. 



4 ACCOUNTING FOR MEASUREMENT 
ERRORS 

Up to this point we did not include measurement errors 
arising from the cluster individual measurements. This is 
however a key issue in the interpretation of the CR-HR di- 
agrams. In this section, we detail our procedure for model- 
ing the measurement errors in the synthetic distributions. 
We first describe how count-rate measurement errors im- 
pact the CR-HR diagrams. In a second step, we also esti- 
mate what would be the uncertainties in the cluster mass 
estimates (based on a M-T proxy) for exactly the same set 
of XMM observations. This step is intended to allow us to 
eventually compare the efficiency of the CR-HR method to 
the traditional method based on cluster masses, in the ideal 
case where redshifts are available. 

4.1 Including measurement errors in the CR-HR 
and z-CR-HR diagrams 

CR-HR diagrams involve three measurements for each de- 
tected cluster: the wide-band ([0.5-2] keV) count-rate CRtot, 
and two narrow-band ([1-2] keV and [0.5-1] keV) measure- 
ments, CRi and CR 2 such that HR = CR1/CR2. 

Errors on these measurements mostly come from Pois- 
son fluctuations in the signal, from the background level 
hampering flux measurement up to large projected radii, 
and from the lack of spherical symmetry amplified by PSF 
distortion effects. Errors on a measured quantity X know- 
ing the true underlying value X are expressed through a 



distribution P(X\X). For the purpose of this demonstrative 
paper, we assume a gaussian error model for CRtot, CRi 
and CR2, without bias and having a non-constant scatter of 
the form: 

/T cxpA -Vy CR 

*oa,® = «HiinJ (m^rJ (4) 

For CRtot, CRi and CR2 we assume CRo = 0.03 cts/s 
and (To = 0.003 cts/s (i.e. a 10% relative error in count-rate 
measurement for 300 collected photons). This simple model 
allows to account for the dependence of measurement errors 
on the number of photons as oc x/iV. We checked its validity 
for measurements of CI clusters i n the 10 ks deep XMM-LSS 
field presented in iPacaud et al.l (|2007l ). see Fig. H3 We note 
that including lower flux systems in the diagram would im- 
ply a more precise model describing the increased influence 
of background on these errors. Errors on HR = CR1/CR2 
are estimated by simulating (200) 3 realistic cluster spec- 
tra on a fine (z, Lx , T) grid, then computing their (true) 
CRtot, 1,2 and simulating CR to t,i,2 following the gaussian er- 
ror model presented above. Then, at fixed (CR, HR) values 
we compute the standard deviation <t cr HR |gg from this 
set of simulated values. Figure [7] shows how o" CH |gg and 
°hr|cr hr i m P ac t the predicted CR-HR distribution of clus- 
ters in the sample. As expected from propagating the errors, 
the relative uncertainty on HR is larger than that on CR and 
it increases as the number of collected photons is lower. 

Using this error model, we "blur" the expected 
dn/dCR/dHR distribution, in the same way as a varying 
PSF would affect an image. Practically, this is done by di- 
viding the initial diagram into fine bins, then redistributing 
the information in each bin into its neighbors using a bi- 
dimensional gaussian distribution with scatters cr^gg and 

CT HR|CR,HR- 

In the case where individual cluster redshifts are 
available, a similar procedure is applied to the three- 
dimensional dn/dz/dCK/dHR distribution. As count-rates 
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Figure 4. Dependence of the CR-HR diagram on the six free model parameters. The fiducial model (Table [TJ is represented by the 
black-dashed contours and predicts 570 clusters over 100 deg 2 (10 ks XMM exposure, CI cluster selection). The red-solid contours 
represent the model obtained when one parameter at a time is varied. The corresponding values for each parameter are: f2 m = 0.30 
(+20%), cr 8 = 0.94 (+20%), w = -0.6 (-40 %), 7 Z , M T = -1, 7z,LT = 1 and x c ,o = 0.2 (+100 %). Contour levels stand for the number of 
clusters enclosed by each curve, labelled by steps of 200. The differential quantity dn/dCR/dHR is constant along each contour. 



measurements are independent on redshift precision, mea- 
surement errors consist of two independent components. Er- 
rors on CR and HR are applied on each redshift slice in the 
same way as for the two-dimensional dn/dCR/dHR distri- 
bution. 



Finally, redshift measurement errors are accounted for 
by narrowing or enlarging redshift bins when integrating the 
density in cubic cells (see sect. [5]). 



4.2 Comparison exercise: errors on the estimated 
cluster masses 

A traditional method for the cosmological handling of X-ray 
cluster samples is to compute a mass proxy for each cluster 
and, subsequently, analyze the resulting redshift and mass 
distribution. Since in Sec. [5] we will compare the CR-HR 
method with the traditional approach, we need to model 
the cluster mass accuracy that is obtainable with exactly 
the same X-ray information. Appendix [A] reviews our as- 
sumed procedure for deriving dn/dM/dz, the mass proxy 
being the cluster X-ray temperature. Apart from the intrin- 
sic scatter in the M-T relation, errors on the mass determi- 
nation mainly arise from temperature measurement errors. 
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Figure 5. The CR-HR diagram resolved in redshift bins for the fiducial model. Each contour from the innermost to the outermost 
encloses 10, 20, 30,... clusters detected over 100 deg 2 (10 ks XMM exposure, CI cluster selection). Measurement errors are neglected 
for this figure. Practically, we compute a 3-dimensional differential density dn/dz/dCR/dHR which is then integrated in defined bins 
(AzACRAHR) to provide an histogram of the clusters in the observable space. This representation can be built straight from the 
observable data, without any assumption on cosmological or scaling-laws parameters, and captures the key features of the sample (see 
text). 
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Figure 6. Count-rate measurement errors as a function of the 
measured count- rate, for the 32 CI c lusters detected in the 10 ks 
XMM-LSS field jPacaud et alj|2007l) . The plain black line shows 
the model adopted in Eq. [4] for a 10 ks survey and the dashed line 
is for 40 ks exposure time. 



Considering the parameter set from Table [T] we derive at 
each (z,M 2 ooc) the expected number of photons N phot col- 
lected in the [0.5-10] keV band with a given exposure time 



1.6 




0.01 0.10 1.00 

Count-Rote [0.5-2] keV (cts/s) 

Figure 7. Effect of count-rate measurement errors on the pre- 
dicted CR-HR distribution. Shaded contours: Predicted cluster 
distribution without measurement errors for 10 ks XMM expo- 
sure, CI cluster selection. Unshaded contours: Same distribu- 
tion after taking into account measurement errors on both the 
[0.5-2] keV count-rate and the hardness ratio. Contours enclose 
respectively 30, 60 and 90% of the total expected number of clus- 
ters. 
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Figure 8. Assumed errors on the measured mass lnA/200c as a 
function of redshift and AftoOc- For each cluster we suppose that 
A/200c is obtained by converting the X-ray temperature estimate 
through the M — T scaling law evolved at the cluster redshift. Red 
dashed lines: lines of constant measurement errors on ln(Af200c) 
in the M200C-2 plane, Black dotted lines: net number of photons 
collected in the [0.5-10] keV band used for the spectral fit. Shaded 
contours: fiducial distribution of detected clusters (10 ks XMM 
exposure, CI selection). Each contour encloses respectively 30, 
60 and 90% of the total expected number of clusters. 



T CX p- We compute: 



AI11M21 



AlnT 




fiV P hot\ 


-1/2 AT 


V 400 J 


~T~ 



(5) 



where AT/T|jv=4oo is the relative temperature error that 
would be obtained with a 40 photons spectrum . This quan- 
tity is taken from Fig. Al of IWillis et ail l|2005l ). considering 
their error bars only. 

Figure [8] illustrates our projected errors on mass mea- 
surements for a lOks XMM observation. At high redshift, the 
main source of uncertainty comes from the number of col- 
lected photons, and the more massive the cluster, the better 
the mass measurement. At lower redshifts the relative mea- 
surement error is almost independent on the mass for the 
clusters being studied. This is because the temperature of 
massive, hot clusters is mo re difficult to dete rmine as they 
lack emission features (e.g. IWillis et~alll2005l ). 

Assuming redshifts are poorly determined - i.e. if only 
photometric redshifts are available - mass measurements are 
further degraded. First, for a pure bremsstrahlung spectrum 
with T(l + z) ^constant at very first order, the spectral 
fit yields a temperature estimate with a relative dispersion 
AT/T ~ Az/(1 + z). Second, the conversion from temper- 
ature to mass depends on redshift through the (1 + z) 7z MT 
factor in equation [2] and a poor knowledge of z impacts the 
mass estimate. We neglect the latter source of uncertainity 
as our fiducial model is computed at 7 z ,mt = 0. The former 
is added in quadrature to the statistical error described in 
eq. [5j, and we take Az/(1 + z) = 0.07 when considering 
photometric redshifts. 



5 THE FISHER ANALYSIS 

In this section, we evaluate the level of performance of 
the method based on the knowledge of dn/dCR/dHR and 
dn/dz/dCH/dRK. We quantify this performance in terms of 
constraints on cosmological parameters and on scaling laws 
related parameters. We describe the Fisher formalism used 
in this analysis and its results. In the ideal case where red- 
shifts are available for each cluster, we compare the efficiency 
of the CR-HR method with the traditional approaches using 
dn/dz and dn/dM/dz. 



5.1 Fisher formalism 

The principle of Fisher matrices applied to 
cosmological fo recasts is thoroughly discussed 



(e.g. iTegmark. Tavlor fc Heavens! Il997l : 



lEisenstein. Hu. fc Tegmarkl 1 19991 : Irleavens! 120091 ). Here 
we briefly recall the approach and show how we applied it 
to evaluate our method. 

Given a set of measured observables {Di,...,D n } as- 
sumed to be uncorrelated, a parametric analysis aims at 
constraining a set of parameters {9i, Op] under a phys- 
ical model M. Defining the likelihood L = P(Di\6^,M) 
and assuming a prior distribution P(# M |.M), the posterior 
P(9fj,\Di,M) oc C x P(9fj,\M) contains all the information 
needed to derive confidence intervals on the If we de- 
note by Oi(^f) t ne observable predicted by the model and 
assuming Poisson distribution in each bin i, the likelihood 
reads: 



In£ = 5ZtoJW(A|Oi 

i 
i 

Defining the Fisher matrix as: 
<9 2 ln£ 



one obtains under those assumptions: 

1 dOidOj 



(6) 



(7) 



(8) 



(9) 



Marginalized parameter uncertainties as well as their 
mutual correlations are encoded in the covariance matrix 
C^v — Fau ■ F° r instance, the 1-a marginalized error on 
parameter # M is given by \J C^. External gaussian priors 
on parameters can be included by simply adding together 
Fisher matrices. Particularly, if Q\ has a prior <7i, the re- 
sulting Fisher matrix is obtained by adding l/crj to the F\\ 
term of the original matrix. 

We insist on the fact that Fisher matrices only provide 
the best constraints attainable by the experiment and ne- 
glect all terms above linear order. The derived constraints 
must thus be seen as indicative for e.g. a comparison of 
two distinct methods. Moreover, a Fisher analysis is valid 
around a given model and all constraints derived from the 
matrix inversion depend on the assumed model. In our case, 
all methods are compared using the fiducial model presented 
in Table [T] Finally, we note that the derivation of eq.|9] pre- 
sented above is valid only if D; does not depend on {# M }. 
In some cases the computation of the data set Di requires 
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the knowledge of some parameters among {6^} (e.g. f2 m is 
needed to compute cosmological distances entering the con- 
version from flux to luminosity). This problem can be par- 
tially overcome by predefining a set {(9 rc f} (it can be the 
same as {^fiducial}), deriving Di with this reference set and 
then correcting Oi(6 ll ) so as to compare both values in the 
same reference space. This is typically the case when deriv- 
ing constraints from the mass distribution of clusters, since 
the mass derivation relies on several key parameters of the 
analysis (see App. . 

Our predicted observable is built from one of the pre- 
dicted densities dn/a!CR/dHR, dn/dz/dCK/dRK, dn/dz 
and dn/dz/dM. Measurement errors are applied following 
the procedure described in Sect. [4] A binning scheme is then 
defined and held fixed, and the Oi are defined as the cell- 
integrated densities. Binning grids are chosen so that the 
bin size at each point is approximately as large as the l-cr 
error size at the considered point. In such a way, correla- 
tions between bins are minimized (we do not consider the 
effect of sample or cosmic variance here, see App. [D] for a 
discussion). We consider two redshift binning : Az — 0.1 and 
Az = 0.03. Observable ranges should span the entire cluster 
population and we choose: z G [0.05, 1.8], CR £ [0.005,3.5] 
cts/s, HR S [0.1,1.55] and Af 2 ooc € [10 13 , 3.10 15 ] ft _1 M . 
Figure [9] shows a typical example of a CR-HR integrated 
density. Comparing it to Figure [4] a substantial amount of 
information has been lost by including measurement errors 
and binning the distribution, but the main characteristics 
of the distribution are still present, in particular its normal- 
ization which is the total expected number of clusters in the 
sample. 

Derivatives of the predicted observables with respect to 
model parameters are evaluated using the five-point stencil 
approximation: 

dO _ 2 O(^ + ^)-0(g,-^) 
96» M ~ 3 89 „ 

12<W M 1 ' 

with steps 5% of the fiducial value for non-zero fiducial 
parameters and 0.05 for the other parameters (7z,mt and 
7z,lt). 

5.2 Results 

In the following, one assumes a 100 sq. deg. survey uniformly 
covered by a 10 ks XMM integration, thus leading to a selec- 
tion function given in Fig. [21 The fiducial model is ACDM 
with parameters given in Table [T] Thus we always consider 
a sample of 570 clusters (see Sect. [3]). 

The analysis involves 15 varying parameters: 

{f2 m , Qb, 0"s, Wo, n s , h}, {QMT, CmT, 7z,MT, ®lnT\M~\ ? 

{olt, Clt, 7z,lt, cr inJj | T } and x c fi- Gaussian priors are 
applied following Table [T] These priors are uncorrelated, 
unless for fib, n s and h for which correlations from WMAP-5 
are taken into account Q. We highlight constraints obtained 
on {n m , o"8, too, 7z,mt, 7z,lt, s c ,o} after marginalization over 
the 9 remaining parameters. 

1 As computed from the the Monte Carlo Markov Chains avail- 
able at http://lambda.gsfc.nasa.gov/ 
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Figure 9. Predicted integrated density dn/dCH/dRH illustrating 
the binning scheme applied in our Fisher analysis (see sect. [5J). 
The original distribution is first "blurred" according to the mea- 
surement errors and then binned into cells. Each bin is approxi- 
mately as large as the measurement error. The model parameters 
are those from Table [l] Shading from white to black represent 
the expected number of clusters in each cell (white = 0, black = 
10 or more clusters). The total number of clusters is 570 over 100 
deg 2 . The error model is defined by equation^ with T CX p = 10 ks 



5. 2. 1 Effect of measurement errors 

In a first step, we study how the precision on count-rate mea- 
surements impacts the constraints on model parameters. We 
consider two situations : i) Count-rates are measured on the 
10 ks survey data and ii) Improved accuracy is provided by 
a subsequent 40 ks X-ray follow-up on each detected cluster, 
(for both cases we have identical selection functions, i.e. the 
same cluster sample). Figure [lOl shows the results for these 
two cases when redshifts are not available, thus using CR- 
HR diagrams only. Constraints on fi m , erg, 7 z ,mt and x Ct o 
show little improvement when dividing measurement errors 
by a factor of two (i.e. going from 10 ks to 40 ks observa- 
tions). On the other hand, constraints on wo and 7z,lt are 
divided by a factor ~ 3. This is a consequence of the defor- 
mation they imprint on the dn/dCR/dHR surface, which is 
better captured in the presence of precise measurements. 

Further, in the case where cluster redshifts are available, 
we investigate the impact of the redshift precision on the 
z-CR-HR method. For this purpose we consider a 10 ks sur- 
vey i) without redshifts, ii) with very approximate redshifts 
and iii) with photometric-like redshift accuracy. Practically, 
these configurations are rendered by narrowing the redshift 
bins from Az — oo, 0.1, 0.03 when computing the Fisher 
matrix. Results are displayed in Fig. 1111 As expected, we 
notice an overall improvement of the constraints obtained on 
model parameters with increasing redshift accuracy. Adding 
redshift information substantially improves the precision on 
wo and 7z,lt, by a factor of 5 (resp. 3); the other parame- 
ters also show an improvement. However, refining the red- 
shift bins does not have a strong impact on the results, and 
a Az = 0.1 binning contains almost the full constraining 
power of the method. 

Table [2] summarizes the 1-a marginalized uncertainties 
on the six parameters in the configurations presented above. 
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Obs. 


dn/dCR/dHR 


dn/dz /dCR/dHR 


dn/d 


z 


Depth 


10 ks 


40 ks 




10 ks 


- 




Az 






0.1 


0.03 


0.1 


0.03 


Param. 
















0.09 


0.08 


0.08 


0.08 


0.15 


0.13 


CT8 


0.14 


0.13 


0.12 


0.12 


0.98 


0.74 


WO 


2.2 


0.73 


0.43 


0.42 


9.0(*) 


5.9(*) 


7z,MT 


0.57 


0.46 


0.44 


0.44 


11(*) 


9.0(*) 


7z,LT 


2.3 


0.78 


0.76 


0.75 


41(*) 


28(*) 


%c,0 


0.04 


0.04 


0.03 


0.03 


0.33 


0.18 



10 ks 
40 ks 



Table 2. Marginalized 1-cr constraints on the six cosmological 
and scaling laws-related parameters for an XMM 100 sq. deg. sur- 
vey at 10 ks depth providing a sample of 570 clusters. We show 
results for three different observables and different measurement 
errors on CR (count-rate), HR (hardness ratio) and z (redshift). 
The 40 ks indication on the second row refers to the depth of a 
potential X-ray follow-up on individual clusters (10 ks meaning 
no additional follow-up). The dn/dz analysis is independent on 
any X-ray follow-up.. The increase in precision of redshift mea- 
surements is rendered by narrowing the redshift bins Az. (*)thcsc 
parameters can be considered as completely unconstrained by the 
observable in this configuration. 



5.2.2 Comparison with dn/dz 

For comparison purpose, we also quote in Table[2]constraints 
obtained from an analysis that would only involve the red- 
shift distribution of clusters in the sample, i.e. the standard 
dn/dz cluster counts not making use of the spectral infor- 
mation potentially available in the X-ray data: even if those 
properties have an implicit impact on the observed dn/dz 
(through the survey selection function), they do not inter- 
fere with the construction of the redshift histogram of the 
sample. Thus, one expects degeneracies between parameters 
to be important and marginalized uncertainties on individ- 
ual parameters to be large, if not meaningless. We discuss 
in App. [C] the comparison between our implementation of 
the dn/dz analysis and a slightly different method based on 
the simultaneous fit of the redshift histogram and the L-M 
relatio n in several redshift bins as presented in lPierre et al.l 
!l20Tll) . 

Figure [T5] illustrates the comparison between this 
method (based on dn/dz) and our method (dn/dCR/dHR) 
which does not make use of individual cluster redshifts. It 
turns out that for dn/dz, parameters u>o, 7z,mt and 7z,lt 
are totally unconstrained and f2 m , as and x Ct o ellipses are 
considerably widened. Moreover, additional degeneracies be- 
tween parameters arise and also participate in diluting the 
constraints. This is particularly true for both parameters 
7z,mt and 7z,lt whose effects cannot be disentangled by the 
redshift distribution only. As expected the best strategy is 
thus to use all information available in the survey (redshifts 
and X-ray measurements) as they help in breaking degen- 
eracies related to the selection function. 



5.2.3 Comparison with dn/dz/dM 

We now assume that redshifts are available for each clus- 
ter, at a sufficient precision to allow a binning size of 0.03, 
corresponding roughly to photometric redshift precision. We 
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Figure 10. Effect of the count-rate precision on the constraints 
obtained by the CR-HR method. Displayed are the constraints 
on the six unknown parameters from Table [T] as obtained by the 
Fisher analysis. Each ellipse encloses the 68% confidence area of 
the marginalized posterior distribution. This figure shows how 
the uncertainty on each parameter can be tightened by reduc- 
ing measurement errors by a factor of ~ 2 (i.e. with a 40 ks 
XMM follow-up of each cluster). On the diagonal is shown the 
marginalized gaussian posterior distribution normalized so as to 
yield a total probability equals to 1. 
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Figure 11. Effect of redshift precision on the constraints ob- 
tained by the (z-)CR-HR method. Two redshift binnings (accu- 
racies) are considered: Az = 0.1 and Az = 0.03. Red ellipses are 
the same as on fig. llOl and show how those constraints are affected 
by removal of the redshift information. 
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Figure 12. Comparison between constraints obtained from 
a traditional dn/dz analysis and our proposed observable 
dn/dCK/dHR, the latter not involving direct redshift measure- 
ments, contrary to dn/dz. The redshift bin size for dn/dz is 
Az = 0.03 and measurement errors are computed for a 10 ks 
survey for dn/dCR/dHR (see fig. 0. Plotting ranges have been 
widened to ease visualization, as in most of the cases the contour 
corresponding to cin/dCR/dHR, is very small in comparison to the 
corresponding green ellipse. 



compare two ways of analysing the data: either directly us- 
ing the observed quantities (dn/dz /dCR/dHR) or using a 
mass proxy (dn/dz/dM) as described in App.0 Neglecting 
measurement errors, intrinsic scatter in the scaling relations 
and systematics, we expect constraints on model parameters 
to be of the same order of magnitude, as both methods use 
identical datasets and rely on the same underlying quan- 
tity (the "cosmological" mass and redshift distribution of 
haloes). Figure [T3l compares the efficiency of the two meth- 
ods, taking into account measurement errors as presented 
in Figures [7] and [S] Corresponding marginalized constraints 
are presented in Table [3] An additional column gives the re- 
sults that would be obtained with a mass precision 10%, 
hence the ultimate constraints attainable with the tradi- 
tional dn/dz/dM function. 

Interestingly, the accuracy reached on parameters of 
purely cosmological origin (f2 m , erg and the dark energy pa- 
rameter wq) is comparable for the two methods. The relative 
precision on Q m is about 30%, while as is constrained to 
~ 15% in all considered cases, showing slight improvements 
when reducing measurement errors and narrowing the red- 
shift binning. On the other hand, a substantial gain is ob- 
tained on both parameters 7 z ,mt and 7 z ,lt governing the 
evolution of scaling laws with redshift. Using the same data 
set, a factor ~ 1.5 appears when using the z-CR-HR rep- 
resentation instead of the z — K'hooc distribution. Finally, 
the parameter x C] o governing the scaling of the /3-model 
core radius r c with i?sooc is constrained twice as well with 
dn/dz/dCR/dHR, up to a relative precision of ~ 30%. 
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Depth 


10 ks 
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Table 3. Marginalized l-cr constraints on the six cosmological 
and scaling laws parameters for an XMM 100 sq. deg. survey at 
10 ks depth providing a sample of 570 clusters. We show results for 
two different observables: (A) is based on the three-dimensional 
dn / dz / dCR / dUR diagrams and (B), the traditional method, re- 
lying on the two-dimensional mass and redshift distribution. The 
survey depth is 10 ks for both the detection and the measure- 
ments. In the last column an uniform mass precision of 10% is 
assumed. The narrowing of the redshift bins Az renders the in- 
crease in precision of redshift measurements. 
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Figure 13. Comparison between constraints obtained from 
our proposed z-CR-HR method (blue ellipses) and from a 
dn/dz/ dM200c analysis (orange ellipses) as from Table \3\ The 
mass M20OC is estimated from the X-ray temperature T then con- 
verted via an assumed M — T scaling relation (see Appendix. [A] 
for details.) The redshift binning is such that Az = 0.03. 

As already mentioned in Sec. 15.2.11 switching from a 
redshift accuracy of 0.03 to 0.1 has not strong impact on 
the final derived constraints for both methods. 



6 SUMMARY AND CONCLUSIONS 

This paper discusses the efficiency of a new method based 
on strictly observable quantities (i.e. instrument dependent 
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measurements) to analyse the cosmological content of large 
X-ray cluster samples. Specifically, for each cluster we only 
make use of an X-ray count-rate, CR, and of a hardness ra- 
tio, HR, plus the cluster redshift, if available ; the 2D (3D) 
distribution of the CR-HR-(z) values from the cluster sam- 
ple constitutes the quantity to be analysed. Compared to the 
traditional approaches based on the dn/dz and dn/dz/dM 
statistics, our method follows a purely top-down procedure, 
requiring to derive from an ab initio model the expected 
CR-HR diagram. The method constrains in a self-consistent 
manner the three main ingredients of the model: (1) the 
cosmology, (2) cluster scaling laws and their evolution and 
(3) the selection effects inherent to the survey design. More- 
over, it avoids the tedious intermediate steps involved in the 
derivation of the cluster mass estimates for the traditional 
methods (for instance : spectral fitting, determination of the 
mean cluster temperature and finally, mass estimate via a 
scaling relation). The present study has been done for the 
particular case of a 100 deg2 survey uniformly paved with 
10 ks XMM observations. We discuss below our results and 
the assumptions made throughout this work. 

6.1 Main results 

- The CR-HR method is primary intended for the analysis 
of X-ray cluster surveys for which no information is avail- 
able on individual cluster redshifts (either from optical spec- 
troscopy or the X-ray observations are too shallow to yield 
an X-ray spectroscopic redshift); this is the case in the early 
phase of surveys covering a large fraction of the sky (e.g. 
pointed observations from archival data). In such a case, the 
usual approach is to consider the logN-logS distribution of 
the cluster fluxes, in one or, possibly, in two or more bands. 
However the constraints one ought to put on the cosmology 
and cluster physics are at this stage, rather degenerate. By 
studying the CR-HR diagrams, we make use of all X-ray 
available information (as far as allowed by the statistical 
significance of the HR), that is: not only we are able to say 
" We have so many clusters in this flux range and in these 
detection bands" , but also ascribe an X-ray colour and coun- 
trate to each cluster, which is much more constraining. At 
the same time, our top-down procedure avoids the non uni- 
versal count-rate => flux translation step, which is mandatory 
for studies based on the logN-logS. 

- The power of the method can be qualitatively intuited 
from Fig. as follows: The traditional way of reading this 
figure leads to the trivial conclusion that redshift is degen- 
erate with temperature, when only an X-ray hardness ratio 
is available: lines of iso-hardness ratio are almost vertical. 
Conversely, one can use this property to infer that the HR 
gives a rough indication of the temperature, independently 
of the redshift; in parallel, for a given temperature the CR 
in the [0.5-2] keV band provides the normalisation of the 
spectrum, which depends on the total cluster emissivity (i.e. 
luminosity) and on the distance of the cluster. Assuming a 
standard M-L relation, it is thus possible to roughly infer 
estimates for the mass and the redshift of a cluster knowing 
the CR and the HR. 

- Practically, our method requires to compute a grid of CR- 
HR diagrams to explore all the parameter ranges one is aim- 
ing to constrain (cosmology, cluster physics and evolution) 
and to find which is the one that best fits the observed dia- 



gram. To evaluate the actual power of the method, we per- 
formed a Fisher analysis which first required a realistic mod- 
eling of how measurements errors on CR and HR dilute the 
information contained in the diagrams. In parallel, for com- 
parison, we performed a similar analysis for a study which 
would be based on exactly the same X-ray data, but would 
determine and use the traditional dn/dz and dn/dz/dM dis- 
tributions. In all of this work, because only a few hundreds 
of clusters are available, we assume that the local M-T and 
L-T scaling relations are known ; but this condition can be 
easily relaxed in the case of all-sky surveys (see paper II). 
We parametrize the evolution of the scaling laws by two fac- 
tors (1 + z) 7 - 

- The calculations presented in this article have been per- 
formed in the case of an XMM/EPIC survey. They can be 
easily extended to any X-ray telescope providing compara- 
ble spectral-imaging capabilities, i.e. a spectral resolution of 
the order of 5-10% between 0.5-2 kev. 

- We summarize our results as follows: 



• The CR-HR method (not requiring redshifts) allows a 
competitive and almost readily available analysis of X-ray 
surveys without the need to wait for spectroscopic follow- 
up of collected clusters. It appears much more efficient than 
the dn/dz statistics (requiring redshift) and thus, is a very 
significant improvement over the logN-logS approach. 

• Refining the precision on CR and HR by multiplying 
the X-ray depth by a factor of 4 (without changing the to- 
tal number of clusters in the sample) does not significantly 
impact the determination of as and Q. m , which are mostly 
dependent on cluster counts. But a significant improvement 
is observed for the cluster evolution and, interestingly, on 
the w parameter of the dark energy equation of state. 

• We further investigated the CR-HR method, by assum- 
ing that redshifts are available for all clusters. This allows 
us to add a third dimension to the diagram: for a redshift 
accuracy of Az ~ 0.1, we observe a significant improvement 
especially for the cluster evolution and the dark energy pa- 
rameters. Increasing the accuracy to dz ~ 0.03 does not 
result in a further improvement, because this is below the 
cluster evolution time-scale. 

• We finally compare the CR-HR-z method to the 
dn/dz/dM statistics. Both approaches appear relatively 
equivalent for the cosmological parameters, while again the 
CR-HR-z method better constrains the cluster evolution pa- 
rameters. 

• In all situations, the CR-HR(-z) approach appears to be 
uniquely well suited to constrain the cluster characteristic 
size. 



In conclusion, the CR-HR method appears optimally 
suited to the analysis of the X-ray data in a survey, and 
the CR-HR-z method is significantly more efficient than the 
standard approaches based on the dn/dz and dn/dz/dM 
statistics and simpler to implement. We have attempted to 
make an account as realistic as possible of the various sources 
of uncertainty entering the method. This is however, proba- 
bly not exhaustive and we discuss further in Appendix |Pl a 
number of pending issues in this respect. 
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6.2 Future work 

In a future work we will study the impact of the shape of 
the selection function on the efficiency of our method. We 
will also quantify the effect of the scatter in the scaling rela- 
tions, especially focusing on the conversion from observable 
to mass in the dn/dz/dM analysis. Another point of interest 
will be to consider how flux measurements in fixed apertures 
can help in a better determination of model parameters, and 
in the case of a much deeper survey, how useful it would be 
to introduce a second hardness ratio pertaining to the harder 
part of the cluster spectrum. In the latter case we shall also 
consider introdu cing Yx as a proxy for the observa ble to 
mass conversion (|Kravtsov. Vikhlinin. & Nagaill2006l ). 

A major, practical, advantage of our method is that 
there is no need to derive individual masses of detected clus- 
ters and we want to investigate how our method can be cou- 
pled to future, full-hydro numerical simulations to constrain 
cosmological parameters without requiring the computation 
of the mass function nor assuming specific scaling laws. 
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APPENDIX A: A REALISTIC MODEL FOR 
THE OBSERVED MASS AND REDSHIFT 
DISTRIBUTIONS 

In this section, we describe our modeling of the derivation of 
dn/dz/dM, in quantity solely used for comparison purpose 
in the present article. In order to realistically introduce the 
mass error measurements in the Fisher analysis, we carefully 
reproduce the various steps involved in the mass determina- 
tion. In this way, we are able to model how the mass accu- 
racy depends on the data, namely, on the number of photons 
collected by the instrument. 

Using the same ingredients as presented in Sect.[3l we 
compute the expected dn/dz/dT/dLx distribution of the 
clusters selected passing the CI selection, which is based 
on the total count-rate and apparent extent. We then as- 
sume that the observer is able to measure the temperature 
T from the collected X-ray photons. The accuracy of such 
a temperature measurement depends mainly on the number 
of photons but also on the cluster temperature and requires 
a prior knowledge of the cluster redshift to a relatively good 
precision. The influence of these factors on the correspond- 
ing observable dn/dz/dM is discussed in Sect. [4] 

The cluster mass is finally obtained by converting 
T via a mass-temperature relation with parameters cho- 
sen in advance. The choice of this "reference" parame- 
ter set is necessary to perform the Fisher analysis de- 
scribed in Sect. [5] In the real data analysis, one may 



also choose to recompute masses with the current param- 
eter values at which the likelihood is est imated (see e.g. 
IVikhlinin et al.ll2009al ; iMantz et al.ll2010al V In practice, we 
integrate dn/dz/dT/dLx over Lx, then at each z we con- 
vert T into M200C using eq.[2]and a "reference" parameter set 
equals to the fiducial model (Table[TJ). We finally redistribute 
the result with a constant scatter a — qmt x o lnT \ M ~ 
a inM\T t° account for the intrinsic scatter in the scaling 
law. Thus in all this paper, the quantity M200C refers to the 
mass obtained from the temperature proxy. The shape of the 
dn/dz/dM distribution for the fiducial model is shown on 
Fig. lAll along with its dependence on the various parameters 
of the model (to be compared with Fig. [4|. 



APPENDIX B: ENERGY RANGES FOR CR-HR 
DIAGRAMS 

This Appendix presents the practical considerations that led 
to the particular choice of the energy bands used in the 
article. 

Figure [Bll displays four synthetic XMM cluster spectra 
for typical temperatures and redshifts. The [0.5-2] keV band 
is optimal for cluster de tection, given the telescope response 
and background levels (|Scharlll2002i ). Moreover, the flux to 
count-rate conversion in this energy range weakly depends 
on the temperature for 0.5 < T < 15 keV and < z < 
1. Consequently, the count-rate in [0.5-2] keV reflects the 
overall normalization of the cluster X-ray spectrum and is 
directly related to the cluster bolometric luminosity. 

A rough estimate of the cluster spectral shape is 
the hardness ratio, basically the ratio between two count- 
rates (or flux measure ments) in two energy bands (see e.g. 
iBohringer et all l200ll for ROSAT clusters and a different 
choice of bands). Because of the particular shape of cluster 
spectra, the high particle background above 2 keV and the 
loss of XMM sensitivity, measurement uncertainties in hard 
bands (typically [2-10] keV) are high. Thus our low-count 
clusters (100-1000) are much better characterized at energies 
below 2 keV. Selected bands must be large enough to min- 
imize the sensitivity to emission features and to maximize 
the signal-to-noise ratio, but sufficiently narrow to be sensi- 
tive to changes in the spectral shape. As shown on Fig. IB1I 
[1-2] and [0.5-1] keV appear as good compromises. 

We note that a deeper exposure (or higher senstivity) 
could allow for complementary measurements in the hard 
part of the spectrum. 



APPENDIX C: CONSTRAINING POWER OF 
THE CLUSTER REDSHIFT NUMBER COUNTS 

We presented in iPierre et all |201ll) cosmological forecasts 
based on the redshift distribution of clusters (and their spa- 
tial correlation function) over 50 deg 2 . To account for the 
unknown scaling relations in the sample, we assumed that 
a mass-luminosity relation can be derived in each of the 20 
redshift bins considered. We parametrized the unknown M- 
L normalizations through 20 parameters at (one per redshift 
bin), further marginalized over when extracting cosmological 
constraints. We did not put any prior on the normalization 
of the scaling law nor assumed a functional form for their 
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Fiducial 
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Redshift Redshift Redshift 

Figure Al. Dependence of the redshift-mass distribution of clusters on the six free model parameters. The fiducial model (Table [TJ is 
represented by the black-dash contours and predicts 570 clusters over 100 deg 2 (10 ks XMM exposure, CI selection). The red contours 
represent the model obtained when one parameter at a time is varied. The corresponding values for each parameter are: f2 m = 0.30 
(+20%), cr 8 = 0.94 (+20%), w = -0.6 (-40 %), 7 Z , M T = -1, 7z,LT = 1 and x c ,o = 0.2 (+100 %). Contour levels stand for the number of 
clusters enclosed by each curve and are drawn by steps of 200. The differential quantity dn/dz/dM is constant along each contour. 



evolution, but we implicitely supposed that the cluster mass 
and luminosity can be derived for each individual cluster in 
the sample (directly from X-ray data or from additional, 
multi- wavelength, observations). 

The present dn/dz analysis differs in that local scaling 
laws are supposed to be known at a fairly good precision, 
requiring the call to an external work. This is rendered by 
putting stringent priors on the parameters by which they 
are defined (Table [TJ. Moreover, the evolution of their nor- 
malization with redshift involves two factors on the form 
(1 + z)t°,mt and (i + 2 )7z,lt^ on the other side, only the 
redshift histogram enters the fitting procedure and there is 



no need to compute physical properties of clusters at any 
moment in the analysis. 

Despite the difficulty of matching Fisher forecasts ob- 
tained by different modeling, we performed a comparison 
between both approaches. In the former one, we let f2 m and 
as free in the Fi sher analysis, while priors on the aj's were 
set according to IPierre et al.l (|201ll ). i.e. assuming a mass 
accuracy of AlnM = 0.5 (corresponding to a 10 ks XMM 
exposure, CI cluster selection). In the current analysis, we 
put priors on the normalization of the L-T relation as well 
as on its evolution parameter such that <t(Clt) = 0.02 and 
o"(7z,lt) = 0.5 with a correlation of p(Clt,7z,lt) = —0.9 
in order to mimic the priors on the ai's that would come 
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Channel energy (keV) 

Figure Bl. Synthetic APEC spectra convolved by the XMM 
response, as they would be observed at a very high signal-to-noise 
level. Vertical lines correspond to 0.5, 1 and 2 keV, i.e. the limits of 
the 3 energy bands of interest. Each of these spectra is defined by 
z, T (keV) and Lx (10 44 ergs/cm 2 /s) with the following values: 
Black: (0.3, 1.1, 0.5) ; Green: (0.7, 1.1, 0.5) ; Blue: (0.3, 4.7, 5.1) ; 
Red: (0.7, 4.7, 5.1). 

from the data itself. We found (AOn ~ 0.03, Atj 8 ~ 0.06) 
with the a-method and (Afi m ~ 0.02, Act 8 ~ 0.07) with the 
present method. Keeping in mind the difficulty to account 
for the various parameter degeneracies and the modeling 
differences between both approaches, we conclude from this 
comparison that they are consistent with each other. 



APPENDIX D: INVENTORY OF THE 
SOURCES OF UNCERTAINTY 

In this section we review a number of sources of uncertainty 
and to which extent they were modeled in the present analy- 
sis. Note that most of them are relevant both for the CR-HR 
approach and the traditional methods based on cluster mass 
estimates. 

Dl Cluster scaling laws 

Fitting simultaneously cosmological and scaling laws param- 
eters as considered here, tends to minimize selection biases 
as it fully takes into account the sample selection function 
and provides values of the sc aling laws for a full range of 
cosmological parameters (see iMantz et al.l l2010bl for a re- 
cent application). In particular strong correlations exist be- 
tween cosmological and scaling laws parameters, whatever 
the method used: apart from the well-known f2 m -er g degener- 
acy, figures [10] and [13] show that evolutionnary parameters 
7z,mt and 7 z .lt also correlate strongly with cosmological 
parameters. 

However, our procedure requires prior assumption 
of a model for the M-T and L-T relations. For this 
work, we have chosen two power laws (equations [2] 
and [3]l with constant scatters, m otivated by the observa- 
tion of individual galaxy cluster s (lArnaud fc Evrardl [l99^ . 
lArnaud. Pointecouteau. fe PrattJ l2005l . IPratt et al.1 120091 ) 
and the hydrostatic equilibrium formalism. However, any 
physically motivated modification of the scaling laws could 



in principle be parametrized and studied along with the 
other parameters. In particular, we did not include any evo- 
lution of the scatter in scaling laws, for which observational 
evidence is weak, nor did we introduce intrin sic correla- 
tion b etween luminosity, tempera ture and mass (|Nord et al.l 
120081 . pitm in lMantz et al-lfeoiOal ). For the purpose of com- 
paring different methods, we consider these latter parame- 
ters to have a reduced impact over our results. 

D2 Measurement errors 

Our models of measurement errors intend to include the 
main sources of uncertainty arising in real X-ray cluster 
analyses. Our relative error on the [0.5-2] keV measured 
count-rate amounts to ~ 6 % for a typical 500 counts 
cluster, roughly consis t ent with past and c urrent analyses 
l|Bohringer et all l200ll. iPacaud et all 120071 ) Our assumed 
measurement error on lnM2ooc for a typical 10 14 ft -1 Mq 
cluster at redshift 0.4 (yielding ~ 500 counts with a 10 ks 
XMM exposure) is ~ 0.5, in agreement with iPierre et al.l 
|201ll) . We note that mass measurement errors should in 
principle depend on the assumed cosmology and scaling 
laws. Throughout this analysis, we neglected such variations, 
as we expect this effect to be negligible relatively to the al- 
ready high value of the error. This problem does not affect 
CR-HR diagrams for which no assumption on cosmology is 
needed to derive measurement errors. 

D3 Additional systematics 

D3. 1 Halo mass function uncertainties 

Uncertainties in the predicted cosmological mass function 
is also a source of systematics in real data analyses. In 
particular, they may arise from the different halo-finders 
used by different a uthors to analyse numerical simulations 
l|Knebe et al.ll201ll). a n d amo unt up to ~ 10%. In this work 
we used iTinker et al.l l|2008l ) mass function for A = 200, 
which is calibrated to roughly 5% upon numerical simula- 
tions, provided the cosmological model is close to ACDM. As 
we are comparing constraints from different methods based 
on the same mass function, we expect such uncertainty to 
have negligible impact over our results. We consider this to 
be true for any unaccounted-for systematic error occuring 
before the conversion from halo mass to observables, in par- 
ticular the conversion between different mass definitions. 

D3.2 Profiles and X-ray spatial variations 

Throughout this work we have assumed a very sim- 
ple isothermal, spherically symetric, /3-model with fixed 
P for the cluster X-ray profiles. This assumption en- 
ters the selection function as it is expressed in terms 
of the apparent core radius r c — a; c ,o-^500c It has 
been widely shown that /3-models do not exactly repro- 
duce the actual complexity of X-ray cluster profiles. In 
particular they cannot account for the cool cor e /non- 
cool core discrepancy (see e.g. IPratt fc Arnaudl |2002| . 
ICavaliere. Lapi. fc Fusco-Femiaricll2009l) which can lead to 
selection biases (|Eckert, Molendi. fc Paltanil l201ll ). More 
generally, we neglected spatial variations in the cluster X- 
ray properties and made the somehow strong assumption 
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that a cluster can be described by only three global quanti- 
ties z, T and Lx- We justify this choice by the fact we are 
considering surveys in which the observed X-ray counts per 
cluster collected by the detectors are quite low (between f 00 
and f 000 counts in general) and do not allow for a refined 
morphological analysis. 



D3.3 Total count-rates in wide apertures 

Measured count-rates in the three bands of interest assume 
that the entire cluster profile can be integrated out to a large 
radius independent of the cluster extent, thus neglecting the 
uncertainty due to background misestimation. There are two 
ways of accounting for this systematic: either with large sim- 
ulated samples of realistic c luster observations the n correct- 
ing for the flux loss (see e.g. lBohringer e t al. 2001), or by in- 
dividually fitting a PSF-convolved model o nto the measured 
profil e and integrating it up to large radii ( Barkhouse et al.l 



2006, Pac aud et all 120071 . IVikhlinin et alj 



2009al ). In pres- 



ence of high-quality data, the second option is often pre- 
ferred, although it is model-dependent. We are currently 
investigating how measurements in multiple fixed angular 
apertures can help in improving CR-HR diagrams. Even if 
a model will be needed, we expect it to be parametric and 
"self-calibrated" the same way as we did for x C] o- 



D4 Sample variance 

Throughout this work we neglected the possible correlations 
between neighboring bins and in particular did not take into 
account the sample variance in our analysis. Taking it into 
account would modify th e likelihood expressed in Equ. 
(see e.g. iLima fc Hull2004l l by introducing a covariance ma- 
trix linking the binned observables to each other. The net 
effect of cosmic variance is to lower the constraints on the 
cosmological parameters, but its expression depends on the 
exact shape of the window function. In this work we only 
consider the total surveyed area (100 sq. deg.) without speci- 
fying a survey geometry. We leave this study for future work, 
but we expect the effect of cosmic variance to have the same 
impact over each of our observables as they all derive from 
the same primordial quantity (the distribution in mass and 
redshift) . 

This paper has been typeset from a T^rjX/ JOTgX file prepared 
by the author. 



